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ABSTRACT 


rn  7P  ?!searCh  ?"^0ntrdCt  h°230009  for  the  period  1  September  1972 
to  2,  February  19/3  is  summarized.  A  full  scale  version  of  ar  electro 

data  Pronaneti0UndiHq  Pr°be  -S  described  with  attendent  experimental 

Sia  are°DrespntpHandThCatterinq-ineaSUremen^S  in  limestone  and  dolomite 
cnnsictinnPIf  f Th<:  scatterinq  measurements  are  for  targets 

medium  Measured  J°ints  ar,d  lithologic  contrasts  in  a  dolomite 
S  frequency-dependent  constitutive  parameters  for  the 

imestone  and  dolomite  media  are  given  and  realistic  pulse  propagation 
calculations  using  these  data  are  shown.  The  basic  problem  o?  probe 

described0 1  1S  dlSCUSsed  and  Pr°9ress  toward  solution  of  this  problem 
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calibration  is  discusstd  and  progress  toward  solution  of  this  problem 
described. 


.,1473 


UNCLASSIFIED 

Si-i-u  I  tv  Classification' 


try  words 


Electromagnetic 
Video  pulse 
Hard  rock 
Faul  t 
Contrast 


1 •  TECHNICAL  REPORT  SUMMARY 

i#iii 

i;:s#i:pl  saiis 

^  »  ^Tth"  ISfteS  ?heTd?rb,eu?hf-in^rrogatin9  a  *•*«" 

Sil==piilfe 

?n?e?^rChTht5thnrh"Se"l^?r?“’^ IS *2 

seismic  or 

layering  can  be  "seen”  using  the  system  in  its  other  ™de  ’ 

1)  to  obtai'n'fulTlr^e  ‘he  pros™  during  this  contract  period  are; 
hi,  °P  r  ,  U  scale  field  reasure merits  of  the  type  of  aenlnmVai 
hazards  likely  ro  be  encountered  in  deep  tunnelino  snrh  I, 

131^ 

lit  yCpe  :0n%S  inPg,?hePfr?Prs?ahl0n  M’?  SCatteri" 

has  been  concenfratedTtbe  ,"te",n'-  eT'f°rt 

describedln5 Section"??'  teport  the  purpose  °f  our  research  is 

a  rede  tailed  specific  objectives  of  the  program 

^Ise'^^'ts'Srrffrrei^eSIn^^^^^tSl  “eJ lec‘r8“t?c 

initiation  date  of  the  present  contract,  gi'ven  ?he’ £  locale  " 

’'Sf  dCSCr1bed  '"  Se2l«  »  together 

"  is  also 


1 


signals  into  a  rock  medium*  and  that  the  ocation  of  large  metallic 
objects  on  the  rock  surface  quite  near  the  probe  (the  drilling  machine) 
does  not  affect  the  coupling. 

Section  V  of  the  report  describes  and  illustrates  video  pulse 
propagation  data  taken  in  limestone  and  dolomite  media.  This  section 
also  presents  and  demonstrates  a  technique  for  obtaining  the  frequency- 
dependent  electrical  properties  of  a  medium  from  pulse  propagation  data. 

Successful  electromagnetic  pulse  soundings  of  faults,  joints  and 
lithologic  contrasts  in  a  dolomite  medium  are  shown  in  Section  VI.  It 
is  to  be  noted  that  these  targets  are  clearly  evident  in  the  raw  temporal 
and  spectral  data  -  thus  for  these  targets  a  sophisticated  processing  is 
not  necessary. 

Analytical  pulse  propagation  calculations  are  reported  in  Section 
VII.  These  results  illustrate  the  attenuation  and  dispersion  effects 
of  the  medium  on  video  pulse  signals.  It  is  also  shown  tha"  using  the 
medium  properties  deduced  experimentally  (Section  V)  and  a  state-of-the- 
art  analysis  for  wire  antennas  in  an  arbitrary  medium,  excellent  agree¬ 
ment  between  measured  and  calculated  video  pulse  propagation  is  achieved. 

The  problem  of  probe  calibration,  and  the  necessity  for  this 
calibration,  if  processing  to  obtain  target  impulse  response  waveforms 
is  required,  is  discussed  in  Section  VIII.  It  is  also  shown  here  that 
the  transfer  function  of  the  probe  geometry  in  situ  on  the  rock  surface 
can  be  obtained  experimentally  if  the  rock  medium  offers  a  particular 
geometry.  This  section  also  discusses  analytical  methods  for  modifying 
the  wire  antenna  programs  to  include  an  air-rock  interface. 

Conclusions  regarding  the  electromagnetic  pulse  sounding  system 
are  summarized  in  Section  IX,  and  research  plans  for  the  remaining 
contract  period  are  given  in  Section  X. 


*A1 1  of  the  field  measurements  have  been  made  i r,  limestone  or  dolomite 
quarries.  Dolomite  and  limestone  are  not  igneous-type  rocks  but  their 
electrical  properties  are  quite  similar  to  those  of  an  igneous  rock 
such  as  granite.  To  make  the  measurements  in  igneous  rock  would  have 
required  prohibitive  travel  from  our  Ohio  site.  It  is  felt  however 
that  at  this  stage  the  system  is  fully  ready  for  demonstration  in  a 
true  igneous-rock  tunnel  and  that  the  experimental  results  reported 
herein  fully  justify  this  conclusion. 


In  summary,  significant  progress  in  the  development  of  an 
e  ectromagnetic  pulse  sounding  system  has  been  made.  Measurements  have 
been  madp,  in  the  field,  of  those  geoloqical  targets  defined  to  be  of 
interest  for  a  hazard  detection  device.  It  remains  to  further  measure 
these  same  targets  from  more  difficult  positions  to  explore  the  system' 
full  capabilities  and  to  add  a  final  modification  to  a  theoretical 
analysis  which  will  then  permit  impulse  response-type  processing. 

These  should  be  accomplished  by  the  end  of  the  current  program  At 
that  time,  as  stated  earlier,  it  is  felt  that  the  docurrented  successful 
measurements  reported  herein  fully  justify  a  testing  program  for  the 
electromagnetic  pulse  sounder  in  an  actual  igneous  rock  tunnel. 


1 1 .  PURPOSE 

A  unique  electromagnetic  pulse  sounding  system  using  periodic 
video  pulse  signals  for  the  detection  and  diagnosis  of  hazards  in 
advance  of  rapid  tunneling  operations  is  being  developed.  The 
tunneling  is  assumed  to  be  in  hard  igneous-type  rock  and  the  hazards 
considered  are  planar  and  include  faults,  joints  and  lithologic 
boundaries  of  infinite  extent.  The  research  program  during  this 
contract  period  has  the  following  objectives: 

1)  Obtain  full  scale  electromagnetic  pulse  soundings  of 
faults,  gouge-filled  faults,  joints  and  lithologic 
contrasts  in  a  r^ck  medium  with  electromagnetic  properties 
similar  to  haH  rock.  These  measurements  were  made  in 
Plum  Run  Quarry  near  Peebles,  Ohio  approximately  120  miles 
from  the  ElectroScience  Labordtory.  The  quarry  rock  is 
basically  dolomite  and  the  above  mentioned  geological 
features  exist  in  abundance.  The  quarry  rock  is  a  soft 
rock  medium  but  because  of  the  abundance  of  the  desired 
geological  features  has  been  studied  for  economic  reasons. 
The  essential  difference  between  hard  and  soft  rock  media 
is  the  difference  in  the  conductivity  which  in  essense 
implies  that  the  operating  range  in  hard  rock  will  be 
increased. 

2)  Construct  a  smaller  and  lighter  version  of  the  full 
scale  probe  (antennas)  geometry  which  can  be  easily 
supported  against  a  vertical  rock  face.  The  smaller 
probe  will  greatly  facilitate  certain  of  the  measure¬ 
ments  in  1  and  a  the  same  time  permit  some  experi¬ 
mentation  with  end  terminations  and  synthetic  aperture- 
type  signal  processing. 


3)  Obtain  estimates  of  the  frequency-dependent  relative 
dielectric  constant  and  conductivity  of  the  dolomite 
quarry  rock  via  in  situ  video  pulse  propaqation  measure¬ 
ments.  Certain  measurements  of  this  type  are  useful  in 
interpretation  of  measured  time  and  spectral  tarqet 
signatures . 

4)  Complete  the  theoretical  analysis  and  computer  proarams 
for  adding  an  air-rock  interface  correction  to  a 
previously  developed  computer  program  for  an  arbitrary  - 
wire  antenna  (bare,  insulated  or  partially  insulated) 

in  an  infinite,  lossy,  homogenec us  medium.  This  result 
will  add  a  system  calibration  ard  theoretical  probe 
desiqn  capability  to  an  existinc  capacity  for  realistic 
(frequency-dependent  constitutive  parameters)  pulse 
propagation  calculations. 

5)  Devise  new  signal  processing  schemes  for  electromagnetic 
pulse  sounding  data  to  enhance  the  detection  and  identi¬ 
fication  capabilities  of  the  system.  Basically,  improved 
estimates  of  the  target's  impulse  response  waveform  are 
to  be  obtained.  Given  the  calibration  capabTity  listed 
in  objective  4,  it  will  then  be  possible  to  calculate  the 
transfer  function  of  the  probe  from  measurements  of  a 
known  target  geometry.  Improvec  impulse  response  waveforms 
wi  11  then  be  possible. 


III.  INTRODUCTION 

The  basic  electromagnetic  pulse  sourding  system  has  been  described 
in  previous  reports[l  ,2].  Briefly,  periodic  video  pulse  signals 
with  spectral  content  spanning  the  Hz  to  C Hz  or  Hz  to  MHz  range  are  used 
as  interrogating  signals.  The  interrogating  signal  is  coupled  to  and 
trom  the  medium  by  two  orthogonally  oriented  dipoles  flush  with  the 
medium  surface  with  a  sampling  oscilloscope  acting  as  a  receiver,  "he 
system  operates  in  two  modes;  direct,  with  transmission  and  reception 
on  the  same  dipole  and  orthogonal  with  transmission  on  one  dipole  and 
reception  on  uhe  otner.  In  the  orthogonal  mode,  the  system  is  "blind" 
to  targets  which  are  symmetric  to  an  axis  perpendicular  to  the  dipoles 
and  centered  at  their  feed  points.  Thus  in  the  orthogonal  mode  the 
system  does  nob  see"  the  air-medium  interface  or  similarly  oriented 
layering  wichin  an  isotropic  medium.  The  system  is  automated  i  e 
control  and  processing  of  the  scattered  signal  waveforms  is  via  an’' 
instrumentation  computer.  Extensive  programs  have  been  developed  for 
processing  of  the  temporal  or  spectral  data[2]. 


In  operation  as  a  hazard  detection  device  in  rapid  tunnel  mo 
operations,  the  probe  will  be  mounted  flush  with  the  tunnel  workinq 
faceand  measurements  taken  m  both  operational  modes.  For  tunnels  rf 

processina  win  1°^°  3?  feet)  3  S^nthetic  aperture  !y  e 

L!  fS  " lW  ■  als?  be  used  by  makln9  measurements  at  discrete  * 
points  on  horizontal  and  vertical  diameters  of  the  face. 

that  i!l0the  P^evious  contract[l,2]  it  was  demonstrated  experimentally 
t?  thl  *  IT9^1"9  V-de°  pulse  si9nfl1s  could  be  effectively  coupled 
tu,t'  !hp^rbUrKen  USlng  properly  designed  and  terminated  dfpoles  and 
that  these  probes  were  very  insensitive  to  obstructions  on  the 
a  r-medium  interface.  A  detection  and  identification  capubility  was 

bath  S?  Hrdrd  f°r  3  particu1ar  c1a^s  of  man-made  targets  buried  in 
nnwS  i  and  limestone  media.  These  measurements  were  made  with  low 
power  and  scaled-down  versions  of  the  probe  because  components  for  full 
power  operation  were  not  yet  completed.  ponents  tor  full 

In  this  report,  a  full  scale  version  of  the  electromagnetic 
pulse  sounding  probe  is  described  in  Section  IV.  Tempora,  and  spectral 
-inM  11^trat\n9  the  Performance  of  system  components  are  givenP 
finally  it  is  demonstrated  that  with  proper  design  of  the  dipoles  video 

ri  th^hlVT  bG  Cffectiv0ly  coupled  -  rock  medium  limestone 
near  Jhe  p^be inS6nSltlve  to  1ar9€  batterers  on  the  rock  surface 

Sect' on  V  of  this  report  describes  and  illustrates  video  duIsp 
propagation  data  taken  in  limestone  and  dolomite  media.  Also  in 
Sect^n  V,  frequency-dependent  attenuation  and  constitutive  parameter 
d3J  °bt!lned  the  propagation  measurements  are  given.  Alternative 

Section  vir  dedUCIng  thT  l3tter  Pa,'ameters  are  also  discussed. 

prv-sents  results  of  measurements  on  specific  faults,  joint 

current  contract0  aSt  tar°etS  d°fincd  to  be  of  interest  on  toe 

rock  Jni?e^°n,-VI1’  tb®orePicai  pulse  propagation  calculations  in 
Hip l U  3  3re  flven'  The  effect  of  frequency-dependent  relative 
dielectric  constant  and  conductivity  is  illustrated  Section  VIII 
describes  the  problem  of  probe  calibration  -  a  ne°essary  ^  uU  I 
scattered  field  data  are  to  be  orocessed  to  obtain  thl  LuPl  re- 
sponse  wavefonn  of  a  given  target.  Our  conclusions  with  resrect%o 
SectionmT!^C  puls%sou;dln9  in  rock  media  are  discussed  in 
Section  x*  d  P  anS  f°r  ^  relnaimn9  contract  period  are  given  in 
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IV.  ELECTROMAGNETIC  PULSE  SOUNDING  SYST! M 

sket^'  of.the  dipole  probe  geometry  and  block  diaqram  of  the 
s°u"d’"?  Sf ten  w,th  the  computer  on-line  is  shown  in  Fig?  ™a  and  lb 
respectively  In  Fig.  la,  the  transition  from  coaxial  able  to  300 

S'balanMdTo  ?nnV1?  T?1a’  wl9h  power  b,'oadba"d  bal“"5  (50  ohm 
u.  baTanced  to  200  oh.n  balanced)  specificali  <  built  to  handle  the  two 

do,j;se  generators  used  with  the  system.  The  purpose  of  the  * 

twin  lead  is  to  assure  a  balanced  input  to  the  antenna  by  use  of  a 

forethe°twoddi’nn?r  ^  iS  balanced>  t0  Physically  separate  che  baluns 
reIo2M  h!  P  thUS  P^event^n9  coupling  across  the  baluns  and  to 
remove  the  balun  response  from  the  time  rfgion  of  interest  A  second 

PP™  ee  r  ioU  CU  lyfem  15  S?rn  ,n  Fi"  U-  Thismu^smal  °"d 
probe  in  Fig.  lc  is  fed  in  exactly  the  sane  way,  i.e.  cable  balun 

and  twin-lead  line.  With  the  small  probe,  a  2  inch  layer  of’hairflex 

absorber  is  inserted  between  the  probe  and  the  air-ground  interface 

cObrse ,  different  baluns  and  shorter  lenqths  of  twin  lead  are  us^d 

The  probe  shown  in  Fig.  lc  was  developed  on  another  program  for  use 

Tn  Jh  nV6ry  narr°W  Vlde°  pu1se  generator  (Ikor  generator  -  see  below) 
in  shallow  overburden  measurements.  Its  small  size  proved  convenient 
for  certain  propigation  tests  and  later  for  fault  delineation  and  it 
was  made  available  for  our  use.  However,  design  d„ta  and  tests  Sf  tie 

patent^app^icaUons?6  pr°pHetar>-  "Ub  tbat  b™9™"  aP°^or  pending 

Two  video  pulse  generators  are  used  with  the  pulse  sounding 

Thptnm’-a-?eW ]U  Packard  model  214A  and  an  Ikor  model  R-100. 

The  nominal  output  of  these  generators  is  respectively  a  50  volt  Deak 

2I0V TdBW: ?S  b  repet’F?"  rate  <Hp)  a"d  a  1000  von  eak 

of  the  HP  '  ?fh*  0  Hz  rate  (Ikor).  The  output  pulses 

d;'^Ikor  generators  are  shown  in  Figs.  2a  and  3a  respectively 

b  themself  ?htte"uation  inser:ed  to  protect  the  sampling  head. 

,th®  ca*e  of  the  Ikor  generato  ’  ^Fig.  3)  the  waveform  in  Fi q  3a 

cab  e fS  ThiseueCtS  °f  the.feed  ' cables  (*00  feet  of  RG-9  coaxial 
t-nnl'  •  V*  necessary  in  order  to  obtain  sufficient  delay  -  the 

external  1  v^ic  the  ,Ikor  generator  (which  cannot  be  triggered 

^anT^^fe^?  £?h  cables  ^hTl  ^  Fl'9S- 

Q°I  So  ic'd  are  S"°Wn’  i-6-’  the  prabe  i>s  reTOve^an^thehsections1  °nS 
that  the  ad  are  shorted  together.  The  corresponding  amplitude  (note 
that  the  amplitude  spectra  are  normalized  by  setting  the  larqest 

Fi*qs  1  anqP3'atUde,h°  un!ty^.and  phase  sPectra  for  the  waveforms  in 
Figs.  2  and  3  are  shown  in  Figs.  4  and  5,  respectively.  For  the  results 

s  own,  the  twin  lead  was  lying  on  a  limestone  medium.  The  severe 
attenuation  of  the  high  frequency  content  of  the  Ikor  pulser  (Fi q  5)  is 

Ssed^^t  hioh  f^P  °n9  •en9th  0fhntwin  lead  transmission  line  (309feet) 
used  At  high  frequencies,  roughly  above  100  MHz,  the  300  ohm  twin 

ead  lying  on  a  surface  is  extremely  lossy.  For  example,  if  the  twin 

ea  1s  c'hortened  to  a  total  length  of  6  feet,  the  improved  amplitude 
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Fig.  5.  Spectra  of  Iko>  pulse. 

(a)  amplitude 

(b)  phase 


spectrur'  shown  dashed  in  Fig.  5a  results.  With  the  Ikor  pulse  generator, 
a  physical  separation  of  the  baluns  (at  least  2  feet)  is  necessary  to 
prevent  coupling  across  the  baluns.  It  is  also  desirable,  with  either 
pulser,  to  remove  the  balun  response  from  the  time  region  of  interest. 

A  slight  improvement  in  the  attenuation  shown  in  Fig.  4  could  be  obtained 
by  reducing  the  15  foot  lengths  of  twin  lead.  This  length  however  per¬ 
mits  the  balun  holder  to  be  firmly  supported  by  polypropelene  lines 
which  are  used  across  the  quadrants  of  the  dipoles  to  rigidize  the  large 
probe  in  Fig.  1  when  it  is  oriented  in  a  vertical  position  against  a 
rock  face.  For  this  reason  no  attempt  to  reduce  the  twin  lead  lengths 
was  made.  Much  the  same  comment  applies  to  the  length  of  feed  (coaxial) 
cables.  The  nominal  length  of  the  transmit  and  receive  cables  is  200 
feet.  In  many  field  situations  it  would  be  possible  to  reduce  these 
lengths  to  10  feet  or  less,  thereby  again  decreasing  the  attenuation. 

In  other  situations  however  the  full  length  of  the  cables  is  needed.  It 
is  not  practical  to  continually  alter  the  cable  lengths,  therefore  the 
lengths  were  fixed  in  all  measurements.  *0ver  the  period  when  active 
measurements  were  in  progress  in  the  limestone  quarry,  the  changes  seen 
in  monitoring  response  waveforms  similar  to  Figs.  4  anc  5  were  insufficient 
tc  warrant  recording.  This  is  not  to  say  that  seasonal  vacations  would 
not  occur.  Note  however  that  the  only  component  susceptible  to  changes 
in  the  electrical  properties  of  the  medium  is  the  short  length  of  twin 
lead  -  and  indirectly  the  characteristic  impedance  seen  by  the  baluns. 

This  type  of  measurement  will  be  repeated  and  monitored  during  the  sumner 
measurement  period.  Some  changes  are  anticipated,  particularly  since 
both  pulse  generators  have  been  returned  to  the  manufacturer  for  repairs 
since  our  last  series  of  measurements. 

The  actual  pulse  transmitted  into  the  medium  is,  of  course,  an 
entirely  different  matter.  It  is  not  possible  to  measure  this  wave¬ 
form  and  a  system  calibration  is  required  in  order  to  obtain  it.  This 
point  will  be  discussed  later.  The  two  pulse  generators  whose 
characteristics  are  shown  in  Figs.  2  and  3  provide  different  capabilities. 
The  very  narrow  pulse  in  Fig.  3  affords  a  fine  resolution  for  close-in 
targets  and  yields  a  wide  spectral  coverage  but  the  power  at  any  given 
frequency  is  actually  quite  small  (<20mW/MHz).  The  broader  pulse 
(Fig.  2)  has  much  less  resolution  and  a  narrower  spectral  coverage,  but 
does  yield  significantly  greater  power  within  its  spectrum  (>110mW/MHz). 
Roughly,  the  broad  pulse  has  a  resolution  capability  of  6.75/n  meters  and 
the  narrow  pulse  a  capability  of  0.075/n  meters,  where  n  is  the  effective 
refractive  index  for  the  medium.  The  above  numbers  are  lower  bounds 
which  neglect  any  dispersive  effects  of  the  medium. 

Operation  of  the  sounder  can  be  described  qualitatively  as 
follows:  Periodically,  a  wave  is  launched  on  the  transmit  arm  of  the 
dipole  at  the  feed  point  which  penetrates  into  the  medium.  As  the  wave 
propagates  along  the  dipole  its  fringe  fields  penetrate  deeper  into  the 
medium.  For  shallow  targets,  the  dipole  is  not  currently  being  employed 


as  a  radiating  structure  in  the  normal  sense  since  it  is  the  near-fields 
which  are  being  used.  For  the  deaper  targets,  the  dipole  is  acting 
as  a  radiator.  Reflected  signals  from  the  medium  via  the  same  or 
orthogonal  dipole  are  observed  on  a  sampling  oscilloscope  acting  as 
a  receiver.  he  measurements  are  automated,  i.e.,  the  computer  controls 
the  oscilloscope  sweep  and  records  the  response  waveform. 

A.  Dipole  Design 

The  design  of  the  dipole  structure  for  soundings  of  =>  specific 
medium  has  two  primary  objectives;  to  couple  as  much  of  the  available 
electromagnetic  energy  into  the  medium  as  possible,  and  to  provide  the 
clearest  "time  window"  possible  by  minimizing  the  reflections  from  the 
balun,  the  dipole  feed  point,  the  ends  of  the  dipole,  the  medium  surface 
and  obstructions  on  or  above  the  surface.  The  objective  is  to  isolate 
in  time  the  scattered  pulse  associated  with  a  subsurface  target  from 
the  above  noted  clutter  mechanisms. 

The  time  domain  design  of  antennas  requires  somewhat  different 
considerations  than  conventional  frequency  domain  approaches.  It  is 
convenient  to  view  the  structure  (include  feed  cables,  baluns  etc.)  as 
a  transmission  line  with  changing  characteristic  impedances  along  the 
line.  Each  discontinuity  in  characteristic  impedance  is  examined 
individually,  and  minimized,  without  regard  to  any  other  point  along 
the  line.  If  an  experimental  approach  is  taken,  as  in  this  case,  then 
a  very  long  pulse  (iisecs)  rather  than  a  video  pulse  is  used,  since 
reflections  from  the  leading  and  trailing  edges  of  the  pulse  tend  to 
become  confusing.  With  the  long  pulse,  the  trailing  edge  of  the  pulse 
is  far  beyond  the  time  region  of  interest.  Thus  the  antenna  structure 
is  interrogated  in  the  direct  reflection  mode  by  a  step  function  with 
a  finite  rise  time.  The  rise  time  is  determined  by  the  resolution 
required,  i  e.,  the  size  of  the  antenna  structure  and  spacing  be twee $ 
discontinuities.  A  typical  measurement  is  shown  in  Fig.  6,  where  the 
interrogating  step  has  a  15  ns  rise  time  and  the  dipole  is  that  shown 
in  Fig.  la.  Dashed  lines  at  plus  and  minus  600  mV  correspond  to  open 
and  short  circuits,  respectively  at  the  input  port  of  the  balun,  i.e., 
at  the  end  of  the  coaxial  cable.  Reflections  at  the  balun,  the  feed 
point  and  tne  end  of  the  dipole  are  clearly  seen.  Fig.  6  also  illus¬ 
trates  the  control  of  the  end  reflections  affected  by  varying  amounts 
of  conducting  sheet  (foil)  across  the  ends  of  the  dipole  arms.  In  Fig. 
7,  control  on  the  feed  point  •'eflection  effected  by  changing  the  shape 
of  the  dipole  arm  (and  therefore  the  characteristic  impedance  at  the 
feed  point)  is  shown.  In  chi s  case  the  medium  is  the  overburden  and 
the  linear  antenna  element'  .re  plastic  coated  AWG12  solid  copper  wire. 
Again  the  incident  pulse  has  a  15  ns  rise  time  which  is  evident  from 
the  open  and  short  circuit  reflections.  For  the  measurements  shown  in 
Fig.  7,  the  section  of  twin  lead  was  removed,  i.e.,  the  baluns  were 
connected  directly  to  the  feed  terminals.  It  is  evident  in  Fig.  7  that 
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Note  that  from  Fig.  7  one  can,  by  estimating  the  reflection 
coefficient  at  the  feed  point,  obtain  the  characteri stic  impedance 
of  the  various  dipole  arm  configurations  at  the  feed  point  The 
characteristic  impedances  are  250,  262  and  400  ohms  respectively  for 
: a^ns1n9l0  Vlre  configurations,  for  the  antenna  of  Fig.  la 
with  the  300  ohm  twin  lead  in  place  it  is  necessary  to  obtain  short  and 
open  circuit  waveforms  at  both  the  balun  output  and  dipole  feed  point 
From  these  results,  the  characteristic  impedance  of  the  300  ohm  twin 
lead  lying  on  a  limestone  medium  is  268  ohms  and  that  of  the  dipole  on 
limestone  is  161  ohms.  Some  improvement  could  be  obtained  by  designing 
components  specifically  for  this  task  rather  than  using  "on  the  shelf" 
units.  Results  for  a  dolomite  medium  were  approximately  the  same, 
indicating  that  the  electromagnetic  properties  of  lirrestone  and  dolomite 
are  quite  similar. 


the 


Note  that  in  Figs.  6  and  7  (using  the  direct  reflection  mode) 
input  signal  was  quite  low  to  protect  the  sampling  head  which  has  a 
^  volt  limitation.  A  voltage  limiter  designed  to  protect  the 
oscilloscope  sampling  head  from  harmful  high  voltage  pulses  has  now  been 
designed  ani  built  which  permits  direct  reflection  mode  operation  of  the 
underground  probe  system.  A  block  diagram  of  the  layout  is  shown  in 
'. ig‘  8a‘  direct  reflection  mode  of  operation,  a  single  antenna 

is  used  as  both  transmitter  and  receiver,  so  that  the  samplinq  head 
must  obviously  be  exposed  to  the  high  voltage  input  pu’se.  A  pair  of 
step  recovery  diodes  (SRD's),  which  are  used  as  charge  controlled 
switches  provide  the  limiting  action  necessary  to  protect  the  sampling 
head.  When  charge  is  inserted  into  an  SRF  by  forward  bias,  the  diode 
appears  as  a  very  low  impedance  (<lp). 
removed,  the  diode  continues  to  appear 
the  stored  charge  is  removed,  at  which 
a  high  impedance.* 


When  this  charge  is  being 
as  a  low  impedance  until  all 
point  it  quickly  switches  to 


In  this  application,  the  oncoming  hi qh  voltage  pulse  provides 
the  forward  bias  to  an  SRD  shunted  across  a  small  section  of  50  ohm 
microstrip  transmission  line  (Fig.  8b).  When  the  voltage  of  the  input 
pulse  becomes  greater  than  roughly  800  mV,  the  diode  switches  to  a 
low  impedance,  causing  the  limiter  to  appear  as  a  short  circuit, 
t  ereby  reflecting  the  incident  pulse  back  down  the  line.  After  the 
trailing  edge  of  the  pulse  has  been  reflected,  the  diode  remains  a  low 
impedance  until  all  of  the  stored  charge  is  removed,  at  which  point  it 
switches  back  to  a  high  impedance,  which  is  unnoticable  in  parallel  with 
the  50  transmission  line.  Any  signal  that  is  less  than  800  mV  in 
magnitude  as  most  target  returns  happen  to  be,  passes  through  the 
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(c)  LIMITER  (SIDE  VIEW) 


Limiter  for  sampling  head  protection 

(a)  block  diagram  with  limiter 

(b)  detail  of  limiter 

(c)  detail  of  limiter 


limiter  undistorted.  Actually,  only  one1  SRI.)  is  required  to  protect 
the  sampling  head  from  the  input  pulse,  however,  for  added  protection 
from  reflections  from  accidental  short  circuits  or  large  magnitude 
target  returns,  a  second  diode  was  inserted. 

ligure  8a  shows  the  direct  reflection  mode  system  with  the 
limiter.  The  length  e,-|  can  be  made  as  short  as  possible,  but  zo  and 
0.3  must  be  long  enough  (have  sufficient  time  delay)  so  that  reflections 
from  the  tee  junction  do  not  clutter  the  time  window  of  interest  The 
stripline  limiter  simply  consists  of  two  parallel  plates  with  a  center 
conducter  sandwiched  in  a  dielectric  between  them  (Fig.  8b).  The 
SRD's  are  shunted  between  the  center  conducter  and  the  parallel  plates 
as  shown  in  Fig.  8c.  The  diodes  are  Hewl i tt-Packard  Type  5082-0300 
with  a  maximum  breakdown  voltage  of  100  volts  and  a  maximum  transition 
time  of  360  psec.  These  specifications  make  the  diodes  ideal  for  use 
with  the  50  volt,  45  nsec  HP  pulse  generator. 


The  waveform  shown  in  Fig.  9  clearly  shows  the  relationship  between 
the  time  delays  involved  and  the  lengths  v2  and  ?.3.  The  operation 
proceeds  as  follows:  A  pulse  is  generated  and  propagates  down  the 
line  «.•]  to  the  tee  junction.  The  impedance  looking  into  the  junction 
is  25:  ,  (two  50s)  cables  in  parallel),  so  a  portion  (p=  -1/3)  of  the 
voltage  is  reflected  back  to  the  pulse  generator  where  it  is  dissipated 
in  the  matched  50  output  impedance  of  the  pulse  generator.  The 
remaining  portion  of  the  incident  pulse  i<  split  equally  between  the 
two  50  cables  ?,2  and  f.3.  Assuming  a  50  volt  input  pulse  (HP)  and 
neglecti ng^attenuat ion  of  the  cables,  a  voltage  reflection  coefficient 

of  p  =  25+50  =  3  at  the  tee  imP1ies  that  about  24  volts  is  being 

applied  to  the  antenna  and  limiter.  The  first  pulse  shown  in  Fig.  9 
is  this  24  volt  input  pulse  limited  to  approximately  800  mV.  The 
second  and  third  pulses  are  the  reflections  from  the  balun  and  from 
the  antenna  and  ground,  respectively.  The  antenna  was  not  located 
near  a  target  at  this  time,  so  there  is  no  target  signal.  The  time 
delay  (2t2)  from  the  input  pulse  to  the  balun  reflection  corresponds 
to  twice  the  time  J^lay  of  line  z2.  The  remaining  pulses  in  Fig.  9 
are  multiple  reflections.  The  second  balun  reflection  and  the  second 
antenna  reflection  result  from  the  first  balun  and  antenna  reflections 
having  been  reflected  back  down  z2  to  the  antenna  from  the  mismatch 
at  the  tee  junction  and  then  back  to  the  tee  and  down  (*3)  to  the 
sampling  head.  The  last  pulse  is  a  second  reflection  of  the  first 
pulse  which  had  been  reflected  back  along  line  £3  by  the  short-circuited 
limiter  to  the  tee  junction  which  in  turn  partially  (p=  -1/3)  reflects  a 
pul^e  back  to  the  sampling  head.  This  method  admittedly  can  be  more 
confusing  than  the  orthogonal  mode  because  of  the  multiple  reflections, 
and  slightly  inefficient  because  of  the  signal  power  lost  in  the 
ref  ections  at  the  tee  junction,  but  the  clutter  problem  caused  by  the 
reflections  can  be  overcome  by  proper  choice  of  the  lengths  z2  and  ?3, 
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alternatively,  the  multiple  reflections  could  be  eliminated  by  matching 
the  system  with  respect  to  the  reflected  signal  port.  A  25  ohm  series 
resistor  placed  in  the  £2  port  would  provide  such  a  match.  The  reflec¬ 
tion  losses  are  overshadowed  by  the  advantages  of  working  with  only 
one  antenna  and  the  elimination  of  the  direct  coupling  between  two 
antennas  in  the  orthogonal  mode.  It  is  interesting  to  note  that  a 
limiter  capable  of  handling  both  pulse  generators  is  now  available 
commercially.  This  limiter  has  been  purchased  on  another  program  but 
should  be  available  part  time  for  our  use  during  the  summer  measure¬ 
ment  period.  Thus  it  will  be  possible  to  make  both  direct  and 
orthogonal  mode  measurements  of  various  targets.  The  limiter  described 
here  was  not  completed  in  time  to  be  used  when  earlier  measurements  were 
made . 


It  was  noted  earlier  that  the  probe  configuration  was  very 
insensitive  in  tf.’  orthogonal  mode  to  lar<ie  scatterers  lying  on  the 
surface  when  un  the  overburden.  Similar  observations  can  be  made  with 
regard  to  the  large  probe  on  a  limestone  or  dolomite  surface,  and  this 
is  with  full  power,  with  and  without  a  station  wagon  parked  within  one 
of  the  quadrants  of  the  dipole  arms  for  both  pulse  generators  (see  Fig. 
31).  No  discemable  change  in  the  respon'e  waveform  is  seen.  We  conclude 
therefore  that  the  dipole  geometries  can  be  well  matched  to  rock  or  soil 
media  (without  ground  rods  for  the  case  of  rock),  and  are  insensitive  to 
targets  on  the  surface.  Thus  in  a  tunnel,  the  presence  of  large  equipment 
immediately  behind  the  working  face  should  not  interfer  with  probe 
operation.  One  might  infer  from  this  that  the  sides  of  the  tunnel 
will  also  not  interfere.  This  has  not  been  tested  as  yet,  but  will  be 
during  the  summer  measurement  period. 


V.  PROPAGATION  MEASUREMENTS  AND  CONST FUTIVE 

PARAMETER  ESTIMATES 

In  order  to  interpret  pulse  sounding  data  for  a  particular 
medium  and  to  gate  the  response  waveform  to  interrogate  particular 
depths,  it  is  necessary  to  obtain  the  effective  refractive  index  for 
that  medium  and  for  the  particular  pulse  generator  being  us*d.  Clearly, 
the  effective  refractive  index  for  the  medium  depends  upon  both  the 
spectral  content  of  the  pulse  and  the  propagation  path  length.  Assume 
that  pulse  propagation  measurements  are  made  for  path  lengths  £  and 
J+A£,  where  A£  is  sufficiently  small  that  the  change  in  path  length 
is  the  only  change  in  the  experiment.  Taking  the  Fourier  transform 
of  both  received  pulses  and  normalizing  the  longer  path  by  the  shorter 
one  obtains 
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result*,  for  the  smaller  term  in  the  expression  for  the  complex  propa- 
gatior  factor.  At  the  same  time,  however,  if  at  low  frequencies  the 
relative  dielectric  constant  has  little  or  no  effect  on  the  pulse 
propagation  then  it  is  really  of  no  interest  for  our  purpose.  The  same 
is  true  of  very  low  conductivities,  where  propagation  path  lengths  would 
be  quite  long  before  the  conductivity  effects  became  noticeable.  There¬ 
fore,  the  number  of  practical  limitations  on  obtaining  relative  dielectric 
constant  and  conductivity  data  from  pulse  propagation  measurements  would 
not  appear  to  preclude  use  of  the  method  for  our  purpose. 

Fig.  10  shows  a  sketch  of  a  tunnel  geometry  in  a  limestone  quarry* 
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Fig.  10.  Sketch  of  tunnel  geometry 


larble  Cliff  Quarries,  3135  Trabue  Road,  Columbus,  Ohio 


where  propagation  measurements  were  made  by  transmitting  from  one  of 
the  large  dipoles  of  Fig.  la  on  the  ledge  to  a  small  dipole  (Fiq.  lc) 
held  against  the  roof  of  the  tunnel.  In  Fig.  11a,  the  received  pulses 
for  20  and  30  foot  transmission  paths  are  shown  for  the  45  ns  pulse  of 
Fig.  2a.  The  amplitude  and  phase  spectra  of  the  received  pulses  are 
shown  in  Figs,  lib  and  c.  Corresponding  measurements  for  the  narrow 
pulses  (Fig.  3)  are  shown  in  Fig.  12.  the  effective  refractive 
indices  obtained  were  3.66  and  3.48  respectively  for  the  20  and  30 
foot  path  lengths  for  the  pulses  in  Fig.  11a  and  3.47  and  3.22  for 
those  in  Fig.  12a.  Note  that  the  effective  index  changes  with  both 
pulse  frequency  content  and  path  length.  It  is  apparent  from  the 
spectra  in  Fig.  lib  that  the  processed  results  at  frequencies  greater 
than  50  MHz  will  be  meaningless,  as  is  also  true  for  the  spectrum  of 
the  original  pulse  given  by  Fig.  4.  Much  of  the  character  of  the 
spectrum  of  Fig.  4  is  preserved  in  the  spectrum  of  Fig.  11.  The  major 
difference  lies  in  the  lobe  centered  about  30  MHz.  The  smaller  antenna 
placed  on  the  rock  media  would  be  resonant  in  this  general  region  of 
the  spectrum.  Thus  it  would  probably  be  a  better  radiator  in  this 
frequency  band.  Similarly,  a  250  MHz  limit  seems  apparent  from  Fig.  12b. 
Frequency-dependent  attenuation  and  phase  data  in  dB/foot  deduced 
via  Eq.  (1)  are  shown  in  Fiq.  13.  Figures  13a  and  b  correspond  to 
Fig.  11  and  Figs.  13c  and  d  to  Fig.  12.  The  propagation  data  in 
Figs.  11  and  12  were  taken  in  the  summer  about  a  week  apart  with 
intervening  rains.  It  is  not  surprising  therefore  that  the  attenuation 
levels  are  slightly  different.  No  attempt  was  made  to  obtain  a 
quantitative  estimate  of  the  moisture  content  of  the  limestone.  The 
frequency-dependent  relative  dielectric  constant  and  conductivity  of 
the  limestone  medium  obtained  from  the  propagation  data  in  Fig.  11 
are  shown  in  Figs.  14a  and  14b  respectively.  The  same  parameters 
from  the  propagation  data  in  Fig.  12  are  shown  in  Figs.  15a  and  15b. 

In  Figs.  14  and  15,  available  data  from  the  li terature[4]  are  shown  as 
crosses.  The  moisture  content  is  known  to  have  a  strong  influence  on 
the  electrical  parameters  and  consequently  the  variation  is  not  sur¬ 
prising. 

For  the  dolomite  medium*,  dipole-dipole  type  propagation  measure¬ 
ments  could  not  be  made  at  convenient  ranges.  For  this  medium,  the 
propagation  measurements  wens  made  using  an  inverted  grounded  monopole 
(18"  diameter  circular  ground  plane  and  1  foot  antenna)  and  the  small 
dipole  described  earl ier  (Fig.  lc)  lying  flush  on  the  rock.  A  sketch 
of  the  geometry  is  shown  in  Fig.  16b.  this  geometry  is  not  optimum  for 
this  study  but  represents  a  simple  expedient  for  generating  the  data 
with  the  accuracy  required  for  these  experiments.  Note  that  for  the 
grounded  monopole-dipole  geometry,  k  in  Eq.  (1)  does  not  have  the  same 
simple  interpretation  as  for  the  dipole-dipole  experiment,  i.e., 
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where  e. r  is  the  relative  dielectric  constant  of  the  medium,  o  is  the 
conductivity  of  the  medium,  w  is  the  angular  frequency  and  e0  is  the 
permittivity  of  free  space.  We  have  processed  measured  data"  for  the 
dolomite  medium  as  if  Eq.  (2)  contained  an  equal  sign,  but  know  that 
this  is  a  very  gross  approximation.  During  the  summer  measurement 
period,  dipole-dipole  propagation  data  will  be  taken  using  the  small 
dipoles  flush  with  adjacent  cliff  faces. 

Development  of  a  good  technique  for  measuring,  in  situ,  the 
constitutive  parameters  of  a  medium  over  a  broad  range  of  fiequencies 
is  a  research  problem  of  some  magnitude  in  its  own  right.  Wai t [5 ]  has 
suggested  a  fixed  frequency  two-loop  method  where  the  loops  are  located 
in  a  borehole.  If  the  rapid  excavation  machirery  uses  a  pilot  bore  of 
sufficient  siz^  -  the  two-loop  method  extended  to  consideration  of 
pulse  signals  via  fast  Fourier  transform  methods  might  offer  an 
attractive  method  for  monitoring  the  constitutive  parameters  of  the 
rock  mtdium.  Wait  has  also  suggested  a  method  for  measuring  the 
constitutive  parameters  of  a  medium  using  two  small  loops  on  the 
surface  of  the  medium[6].  This  technique  cojld  also  be  extended  to 
pulse  sinnals  and  would  appear  to  be  applicable  within  limits  to  a 
tunnel  working  face.  These  techniques  need  to  be  tested  (as  was  noted 
by  Wait).  If  time  and  funds  permit,  a  two  surface  loop  pulse  measure¬ 
ment  will  be  made  in  the  dolomite  medium  during  the  summer  measure¬ 
ment  period.  It  should  be  noted  however  that  basic  application  of  the 
electromagnetic  pulse  sounder  does  not  require  "a  puiani"  knowledge  of 
the  constitutive  parameters  of  the  medium.  The  fundamental  require¬ 
ment  is  an  estimate  of  the  effective  refractive  index  of  the  medium 
for  a  particular  pulse  signal.  With  this  information,  gating  to 
provide  temporal  and  spectral  signatures  at  particular  depths  can  be 
effected.  As  discussed  in  Section  IX,  a  sophisticated  processing 
would  proceed  to  remove  the  attenuation  and  dispersion  effects  of 
the  medium  via  theoretical  calculations  (for  which  the  frequency- 
dependent  constitutive  parameters  would  be  needed)  and  thus  obtain  an 
impulse  response  waveform  for  the  target  in  question.  Note  particularly 
chat  these  are  not  plane  wave  attenuation  and  dispersion  calculations 
nor  are  they  for  infinitesimal  elements. 

The  time  and  frequency  domain  waveforms  of  the  transmitted 
pulses  of  the  HP  and  Ikor  pulse  generators  for  the  monopole-dipole 
configuration  in  dolomite  are  shown  in  Figs.  16  and  17,  respecti vel v. 
Note  that  over  a  path  length  of  24',  the  received  HP  pulse  (Fig.  16)' 
endures  for  well  over  1  ysec  while  the  received  Ikor  pulse  (Fig.  17) 
endures  for  over  2.5  ysec.  The  excessive  dispersion  is  believed  to  be 
due  primarily  to  the  moisture  content  of  the  rock  -  it  had  been  raining 
for  at  least  24  hours  immediately  prior  to  the  recording  of  these 


measurements.  Poor  coupling  to  the  rock  medium  ma.v  also  be  a  factor 
in  that  the  surface  (2"-3")  of  the  rock  was  actually  a  loose  gravel  and 
di rt  mi xture. 

The  time  domain  waveforms  of  fig.  16  (HP)  were  Fourier 
transformed  to  obtain  the  amplitude  and  phase  spectra  shown  in  Figs. 

16b  and  c  which  upon  subtraction  and  normalization  by  path  length  yield 
the  attenuation  in  dB/ft  and  phase  constants  of  the  rock  medium  as  a 
function  of  frequency  (Figs.  18a, b).  This  process  was  repeated  for 
the  Ikor  (1  ysec)  pulse  waveforms  of  Fig.  17a  to  obtain  the  amplitude 
and  phase  spectra  (Figs.  17b, c)  and  attenuation  (Figs.  18c, d).  From 
the  Ikor  data,  it  is  apparent  that  frequencies  above  6  MHz  were  lost 
in  the  noise  over  a  propagation  distance  of  24  feet,  while  from  the 
HP  data,  frequencies  above  20  MHz  were  lost  over  the  same  distance. 

It  is  important  to  note  that  the  attenuation  vs  frequency  data  in 
Fig.  18  mean  absolutely  nothing  above  these  limiting  frequencies. 

Using  the  Eq.  (1)  assumptions,  the  transmission  data  were 
processed  to  obtain  a  frequency  dependent  relative  dielectric  constant 
and  conductivity  for  the  medium  (Fig.  19).  Again  these  data  should 
be  considered  only  over  the  frequency  ranqe  determined  by  the  significant 
spectral  content  of  the  transmitted  pulses.  Due  to  the  excessive 
attenuation  and  dispersion  and  the  fact  of  Eq.  (2),  these  data  are 
admittedly  very  gross  estimates,  but  even  so,  they  are  still  within 
an  order  of  magnitude  of  published  results[4]. 


VI.  SCATTERING  MEASUREMENTS  IN  DOLOMITE 

An  initial  series  of  orthogonal  mode  measurements  have  been  made 
in  a  dolomite  quarry  (Plum  Run  Stone  Di vi c i on )  near  Peebles,  Ohio. 
Unfortunately,  the  quarry  is  approximately  120  miles  from  the  Electro- 
Science  Laboratory  -  thus  measurements  at  this  site  involve  significant 
travel  costs.  This  quarry  is  the  closest  location  known  where  faults, 
joints  and  lithologic  contrasts  exist  without  an  intervening  overburden 
cover.  A  geological  description  of  the  quarry  was  given  in  the  Research 
Plan  dated  7  September,  1972  ,  and  can  also  be  found  in  Reference  [7]. 

To  accommodate  measurements  at  the  Plum  Run  site,  the  electro¬ 
magnetic  pulse  sounding  system  was  modified  to  permit  operation  from 
a  small  1250  watt  generator.  The  modification  consisted  primarily  of 
the  construction  of  a  small  waveform  generator  to  supply  the  ramp  and 
gate  waveforms  necessary  for  remote  site  operation.  The  entire  system 
can  now  be  transported  in  a  station  wagon  or  similar-type  vehicle. 

Thus  the  large  truck  described  and  shown  in  a  previous  report[lj  is 
no  lonqer  required.  This  also  means,  however,  that  the  boom  on  the 
truck  is  not  available.  As  will  be  seen,  the  boom  would  be  particularly 
advantageous  for  certain  of  the  proposed  measurements.  A  second  probe 
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geometry  was  also  made  available.  This  probe  (shown  in  Fia  Id  was 
developed  on  another  program  and  is  designed  for  remote  probinq  in 

^fle°cet^r1nrofthrc1HtiV?1y  Sha11°W  deP?,'s-  Basical  ly^theend’" 
ret (ect1  ons  from  the  dipole  arms  are  controlled  in  this  case  bv  nWinn 

absorber  between  the  probe  and  the  medium  to  attenuate  mSst  of  the  h  ah 

frequency  energy  before  it  caches  the  ends  of  the  dipo™  The sma 

It°wasWusedef19ned  !°r  opJration  wi th  the  narrow  (Ikor)  pulse  of  Fig.  3 

because  nf  thpP  certa 1.1?.  of,  dolomite  quarry  measurements  primarily 

because  of  the  ease  with  which  it  could  be  set  up  and  roved. 

made  0^  th^f^oo^orth^  results  °f  f  ?rthogonal  mode  measurements 
Li  1  lev-formating  ™  quarry’  wJ1ch  1s  also  referred  to  as  the 

generators  ZLa  thp"?  W3S  ^hat  Fig'  lc  and  tho  Ikor  P^se 

fom  !i?h  th!  n  n'  /h  PP  Waveform  1S  a  control  or  no  target  wave- 

contrast  p  Sth  .  Pear  the  Center  of  the  floor  and  no  known 
contrasts  beneath  the  probe  to  a  depth  of  at  least  35  feet  The  next 

tTJ?vVe^TaShf0W,tthe  effects  a5  tl,e  e™be  geometry  is  displaced 
vertlca  dLnTa  »  aUl  5ccunn9.a'0"9  the  floor.  The  fault  has  a 
up  s  de  of  the  S  Approximately  13  feet  and  the  probe  is  on  the 

see  the  ror^elaHal  f?UU  wa?  reasonably  clean.  One  can  clearly 

the  fsi.lf  between  the  position  of  the  probe  with  respect  to 

the  fault  and  the  corresponding  response  waveform  Note  that  this  is 
U  ing  the  unprocessed  time  waveforms.  In  each  oTth  4  wavefom  s  own 

ff  th  1  ;\hthe  ?ireCt  ““P»"9  across  the  feed  oin?s  ’ 

or  tne  dipoles.  Note  that  the  shape  of  this  leadinq  spike  is  a  r^nii,a 

rouqhlvP8  to  tfT  1nT,rl5a-3b-  In  particular,  the  base"  widih  ?s 

to  the  leve?ness"of  th!',  i1'"  “upling  magnitude  is  quite  sensitive 
CO  tne  i eve  ness  of  the  probe.  One  can  anticipate  therefore  variations 

couJlinqmisnofUnn  thr°UghoUt  a  series  of  measurements.  The  direct 
coupling  is  of  no  consequence  since  it  can  always  be  qated  out  but 

does  serve  as  a  convenient  time  reference  -  that  is  it  explicit-!/ 

DT  “n,n  ‘T  theT  afrival  of  the  interrogating  pulse  at  the  feed 
viola  a  Prod®-  Tbe  propagation  measurements  tak  »n  on  this  dav 
n  f  anwL?S  IVe  de,a*  fo:  tbf  Ike  PIJ1  ser  in  dolomUe  of  ll/ 
and  it  h,HWJ  the  measurements  were  taken  a  light  rain  was  falling 
fir„u  h-?hb^n  ralnlng  steadily  for  ct  least  24  hours.  It  is  dif- 
ficult  with  this  configuration  to  estimate  when  the  response  from  the 
fault  actually  started.  In  the  normal  care,  i.e.,  with  the  tar ^t 
below  the  probe  plane,  it  is  clear  that  tie  shortest  path  is  via  feed 
point  -  target-feed  point  with  the  response  via  the  ends  of  the  di Dole 
arms  arriving  somewhat  later.  Consider  first  the  first  positive  ePa 
and  the  second  negative  peak  in  the  fault  waveforms.  For  all  th™ 
measurements  these  remain  almost  constant  at  delays  of  9  and  15  ns 
respectively  from  the  start  of  the  direct  couplinq  except  at  the  5’ 
foot  distance  where  something  else  has  obviously  been  adoed  to  the 
response  Since  these  characteristics  do  not  rove  as  the  oroSe  is 

^rkd’hU  1S  <rontended  that  they  are  due  to  another  target  beneath  he 
ock  whose  nature  is  unknown.  On  the  other  hand  the  second  positive’ 
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peak  m°ves  out  in  time  from  the  direct  couplinc  and  occurs  respectively 

m^  delay  o?  S/V™  ?1rect  C0Upllnq-  With  thp  estimated  " 
th'it  0the  [  pu!se  ln  doloniite  ™d  considering  a  2-way  path 

from  S  t  etr°  Th  f  ***"  TOVed  resPect1v£ly  2-07  and  6.2  feet  away 
of  2  feet  and  5  feet  iLa?hrf  ?  r!a?°nably  we11  W1‘th  the  known  move, rent 
sketch  in  ??«  ?n  °J  the  1ast  2  measurements  in  Fig.  20.  The 

is  actually  the  stS  nf  1°^^  misleadin9-  The  fault  line  shown 
ib  actually  the  start  of  fracturing  in  the  rock  -  the  actual  vertir^i 

Howevelltlslnlikel  vlh  *t  ‘°  the  H<,ht  0f  the  line  shown- 

taiLtll  tk!  .  y.t  il  resPonse  l>e  attributed  to  the  fault 

taree^  1  Theh^ir??ewhad  be!n  m°Ved  within  range  of  another  subsurface 
target.  The  amplitude  spectra  in  Fig.  21  also  clearlv  disnl^v  thp 

(^stedecS™hinfantth"lenT°mfhredut0  the  n°  tarqet  aiml  1  tu<te  spectra 
reasneo  curve  in  all  three).  The  phase  spectra  are  shown  in  Fin  ?? 

5  H  r?  ij r- 

^idse  spectral  data,  it  is  evident  from  the  phase  spectrum  nf  thP 
foot  waveform  that  the  unknown  target  has  affected  the  phase. 


_n j  tJhunab°ye  set  of  medsurements  were  repeated  with  the  same  probe 

nn  se  gPnerator‘  These  response  waveforms  are  shown  in 

9roet  and  „°0etl™;9'  23  that  vhV°  there  •">  differences  In  te" 

vividly  ?n2?t)ates  thT  T5eT  n0t  Slmp1y  Interpreted.  This 

y  i ustrates  the  advantages  of  a  very  narrow  interment  inn 

tarqetln  l  ei  1"  N°te  a,s°  «» ^thSISJtl  17 S  ofTml 

target  in  the  5  foot  waveform  is  also  evident. 


re  not  possible  because  heavy  rains  caused  standing  water  in  the 
quarry  making  the  base  of  the  vertical  reck  faces  inaccessab  e  It 

period6' t0  makG  theSe  measurements  during  the  summer  measurement 


the  IkorCgeneretorhefurtppSSf°fiIhe  initial  fault  nieasurements  using 
uie  iKor  gene.ator,  furtner  fault  measurements  were  conducted  in  ™ 

attempt  to  locate  a  minor  fault  line  in  a  sectionTf ?Ue nSrrl  that 

had  not  been  accurately  mapped.  With  the  indispensable  help  of  the 

75afoot9wi>deflstriD  thp  fau1 1  ]ocation  was  narrowed  down  to  roughly  a 
toot  wide  strip.  Figure  24  shows  a  sketch  of  the  site  with  the 
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no  target  geometries.  Small  probe  and 
Hewlitt  Packard  generator. 


numbered  X's  marking  the  probe  positions.  As  seen  in  the  figure,  two 
swt  eps  across  the  strip  were  made,  using  the  small  absorber  covered 
probe  and  the  Ikor  pulser.  The  seemingly  random  spacing  between 
measurements  was  necessary  because  of  surface  roughness,  as  it  is 
always  desirable  to  keep  the  probe  as  flat  as  possible  on  the  surface 
of  the  rock.  Figure  25  shows  the  time  domain  waveforms  obtained  at 
the  17  probe  positions  in  Fig.  24  and  an  average  of  all  17  waveforms 
which  is  used  as  a  no  target  reference. 

Observe  the  large  response  in  Fig.  25,  positions  4,  5,  6,  and 
15.  ^  The  response  is  delayed  in  time  as  the  probe  is  moved  away  from 
position  4.  This  trend  implies  that  the  'target"  under  observation  is 
very  close  to  position  4.  While  some  evidence  of  this  "target"  is 
apparent  in  positions  2  and  3,  it  is  rather  difficult  to  associate  a 
similar  correlation  between  the  response  and  distance  from  target  as 
in  positions  4,  5,  and  6  and  Fig.  20.  There  were  also  a  number  of 
surface  cracks  in  the  vicinity,  indicated  by  the  dashed  line  in  Fig.  24, 
which  are  probably  responsible  for  all  of  the  low  magnitude  responses 
in  Fig.  25.  The  target  in  question  here  is  also  difficult  to  describe 
from  a  geometrical  viewpoint.  The  fault  line  need  not  be  perpendicular 
to  the  plane  of  the  surface,  in  fact,  it  may  be  tilted  at  such  an  angle 
as  to  invoke  a  large  response  to  the  probe  in  position  4  and  yet  invoke 
only  a  minor  response  to  position  3.  Note  that  in  the  amplitude  spectra 
of  positions  4,  5,  and  6,  shown  in  Figs.  26  the  characteristic  low 
frequency  fault  response  is  observed,  which  is  similar  to  those  in  Fig. 

21.  Position  15  contains  the  maximum  response  obtained  in  the  second 
sweep  across  the  fault.  Again  the  amplitude  spectrum  of  position  15 
displays  the  same  characteristics.  The  phase  spectra  are  shown  in  Fig. 

27.  These  data,  and  previous  knowledge  of  the  general  direction  of  the 
fault  lead  to  the  conclusion  that  the  fault  under  investigation  follows 
a  line  that  passes  near  positions  4  and  15.  The  precise  location  of  the 
fault,  however,  cannot  be  verified.  This  section  of  the  quarry  has 
since  been  worked,  and  the  head  geologist  generally  agreed  with  our 
conclusions.  However,  he  stated  that  so  much  fracturing  occurred  in 
this  area  that  precise  location  was  not  possible. 

Perhaps  the  most  striking  results  obtained  thus  far  is  the 
response  from  a  lithologic  contrast  at  a  depth  of  approximately  45 
feet.  Figure  28  shows  the  geometry  of  the  measurement.  The  large 
(28'  dipole)  probe  with  the  aluminun  sheets  is  situated  on  a  layer  of 
rock  about  45  feet  above  the  quarry  floor,  and  is  operated  in  the 
orthogonal  mode  with  the  50  volt  45  ns  (HP)  pulse  generator.  Both  rock 
formations  are  dolomite,  but  the  upper  layer  is  an  open  porous  reef- 
type  material  while  the  lower  layer  (quarry  floor)  is  a  dense  material. 
Basically  there  is  a  large  change  in  the  specific  gravity  of  the  material, 
but  the  contrast  is  not  sharp.  The  extent  of  this  gradual  transition 
is  approximately  3  to  5  feet.  On  either  side  of  the  probe  were  joints 
running  vertically  from  the  surface  near  the  probe  to  the  quarry" 
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Orthogonal  Mode  Measurement  on  quarry  floor 
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(b)  amplitude  spectrum 

(c)  phase  spectrum 


fn  .Tbe  shown  in  Fig.  28  is  that  between  the  Peebles 

formation  (cliff)  and  the  Lilley  formaticn  (quarry  floor).  A 
geological  description  of  these  formations  is  given  in  Reference  r 7 1 
Data  to  estimate  the  electrical  properties  of  the  Peebles  formation  ' 
(cliff)  have  not  yet  been  obtained.  Note  that  all  of  the  data  shown 
!-Ct10P  “ere  obtained  during  two  visits  to  the  quarry  of  one 

thP  ,0nHe;5h:  £ompare  h.owever  the  direct  coupling  in  Fig.  29  for 
the  cliff  and  that  shown  in  Fig.  30  for  the  large  probe  usina  the  HP 

generator  ‘WT*  f1?01"'  •W1th  the  ,ar9e  P™'*  a"d  the  Hp 
generator,  the  direct  coupling  is  much  less  sensitive  to  minor 

vanatTons  in  eyelness  than  the  small  probe  with  the  Ikor  generator 

Thus  some  significance  can  be  attached  to  the  change  in  direct  couplinq 

level  -  which  basically  confirms  a  change  in  electrical  properties 

Corings  of  both  formations  have  been  taken  by  the  quarryPpe?sonnel ' 

but  in  any  laboratory  scheme  for  measuring  electrical  properties  w° 

now  have  it  would  be  necessary  to  shape  the  samples  to  existing 

equipment  sizes.  The  quarry  geologist  is  naturally  reluctant  to 

release  core  samples  on  this  basis.  Estimates  of  the  electrical 

?™rtl6S  Wl11  ^  obtained  via  Propagation  measurements  during  the 
summer  measurement  period.  There  is  possibly  a  second  lithologic 
contrast  shown  in  the  quarry  floor  measurement  in  Fig.  29.  The  L i  1  lev 

rouqhl v°40  feet'Tpl1  °°th  *  9  shale  ^pe  Pock  a  depth  of  ^ 

roughly  40  feet  below  the  quarry  floor.  The  response  in  Fiq.  30 

occuring  at  approximately  350  ns  roughly  corresponds  to  this  depth 

out  intemal^fiP  I®U  ln  this  response  since  the  care  taken  to  rule 
out  internal  reflections  was  not  as  extensive  as  those  for  the  cliff 

encouraging :  NeV6rthelcSS  the  pa^'al  confirmation  is  Jlstlnctly 

Fl"ally’  Flp'  31  shows  the  orthogonal  mode  waveforms  with  and 
without  3  car  parked  in  one  quadrant  of  the  probe.  Figure  31a  cor¬ 
responds  to  the  Ikor  pulser  and  Fig.  31b  to  the  HP  pulser.  In  both 
cases,  there  is  essentially  no  change  in  the  recorded  waveform  aqain 
confirming  the  fact  that  the  probe  is  insensitive  to  scatters  on  the 
surface,  even  when  used  over  a  rock  medium. 


VII.  PULSE  PROPAGATION  CALCULATIONS 

Computer  programs  have  been  developed  to  completely  analyze 
arbitrary  wire  antennas  or  arrays  of  antennas  embedded  in  an  infinite 
o  ogeneous  medium[8].  Typical  program  output  consists  of  self  im¬ 
pedance,  mutual  impedance,  current  distribution,  near  zone  fields  and 

d  pripTiPat  ernS<  Fini?e  Activity  of  the  wire  and  insulating 
electric  sleeves  can  also  be  taken  into  account.  A  piecewise 

sinusoidal  expansion  is  assumed  for  the  unknown  cur^  t  distribution 
on  the  antenna,  and  GalerkiVs  method  is  used  to  r  -Ce  the  integral 
equation  to  a  system  of  simultaneous  linear  equations.  AlthouglAhe 


Fig.  31.  Orthogonal  measurements  with  and 
without  surface  obstacle. 

(a)  Ikor  generator 

(b)  Hewlett  Packard  generator 


program  nonmally  operates  in  the  frequency  domain,  speed  and  accuracy 
are  sufficient  to  permit  a  Fourier  Transformation  to  the  time  domain. 
There  is  one  program  limitation  that  is  determined  by  the  amount  of 
computer  storage  available.  The  antenna  must  be  divided  into  a  finite 
number  of  segments  which  is  related  to  the  number  of  terms  in  the 
piecewise  continuous  expansion  for  the  current.  In  order  to  obtain 
sufficient  accuracy,  the  length  of  each  segment  should  be  less  than 
A/4.  At  high  frequencies,  therefore,  the  number  of  elements  could 
easily  become  too  large  for  the  computer  to  accommodate.  The  current 
limit  of  the  Datacraft  6024  ElectroScience  Laboratory  computer  is  50 
elements.  Larger  computers  such  as  the  IBM  370  could  accommodate 
perhaps  250  segments.  The  Fourier  Transform  modification  has  been  added 
to  the  main  program  to  allow  conversion  between  the  time  and  frequency 
domains.  A  time  domain  voltage  pulse  signal  is  applied  to  the  term  nals 
of  the  transmit  antenna,  shown  in  Fig.  32.  This  antenna  is  identical 


Fig.  32.  Antenna  structure  used  for  propagation  calculations. 


to  one  of  the  dipoles  in  Fig.  1c.  Note  that  Fig.  32  shows  the  antenna 
divided  into  14  segments  as  needed  for  the  calculation.  The  radius  of 
the  wires  is  3/8"  and  assumed  to  be  perfect  conductors.  The  input 
pulse  is  the  50  volt,  45  nsec  (HP)  pulse  (Figs.  33a  and  2b)  after  it 
has  been  passed  through  the  cables,  baluns,  and  twin  lead.  It  can 
be  seen  from  the  amplitude  spectrum  of  Fiq.  33b  that  the  most  significant 
frequency  content  of  the  pulse  is  below  50  MHz,  (above  50  MHz  the  ampli¬ 
tude  spectra  is  down  40  dB)  so  the  lengths  of  the  segments  in  Fig. 

32b  had  to  be  less  than  z  =  a/4  =  c/4f/F^  =  4. 45//c^  feet,  which  they  are 
f°r  er<2°:  "^he  input  time  domain  pulse  (Tig.  33a)  is  Fourier  transformed 
and  the  first  50  harmonics  of  the  resulting  complex  ;pectrum  (Fig.  33b) 
are  applied  to  the  terminals  of  the  transmitting  antenna.  An  induced 
voltage  is  calculated  at  the  terminals  of  an  identical  receiving 
antenna  a  specified  distance  away  at  each  of  the  fif.y  harmonics.  The 
receiving  antenna  is  terminated  on  a  300n  load  to  agree  with  the  300fi 
twin  lead  used  in  the  pulse  sounding  system.  This  spectrum  of 
received  voltages  is  then  inverse  Fourier  transformed  to  obtain  the 


received  time  domain  waveform  These  calculations  have  only  been 
carried  out  with  the  50  volt,  45  nsec,  HP  pulse,  because  the  high 
frequency  content  of  the  1000  volt,  250  psec  pulse  results  in  such  a 
large  number  of  segments  that  the  present  computer  system  cannot 

accommodate  the  program.  The  size  of  the  large  antenna  with  the 

aluminum  sheets  also  presents  a  similar  problem. 

Figure  34  shows  the  effect  of  increasing  conductivity  on  the 
transmitted  pulse  for  a  relative  dielectric  consta.it  (er)  cf  10  and 
a  separation  distance  of  20  feet.  As  the  conductivity  is  varied  from 

a  =  .001  mho  per  meter  to  a  =  .01  mho  per  meter,  the  attenuation  for 

plane  wave  propagation  increases  according  to  the  relation 


1/2 


from  .03774  to  .16953.  In  dB  per  meter,  this  corresponds  to  roughly 
.3  dB/m  and  1.5  dB/m.  The  changes  in  attenuation  can  best  be  seen  by 
comparing  the  amplitude  spectra  of  Fig.  34b.  Note  that  for  very  low 
frequencies  2  MHz),  the  received  voltage  increases  with  increasing 
conductivity,  but  for  high  frequencies,  the  received  voltage  decreases 
with  increasing  conductivity.  Note  also  from  the  phase  spectra  of 
Fig.  34c,  that  the  phase  which  is  rotating  in  a  counter  clockwise 
direction,  is  further  delayed  by  increasing  conductivity.  The  upper 
waveform,  calculated  with  a  =  .001  mho/m,  is  approaching  the  type  of 
waveform  expected  for  a  hard  rock  medium.  Typical  values  for  the 
conductivity  of  a  hard  rock  such  as  granite  are  less  than  .001  itfio/m, 
so  the  magnitude  of  the  transmitted  pulse  can  be  expected  to  increase. 

The  ringing  associated  with  the  waveforms  for  low  values  of  conductivity 
can  be  decreased  by  adjusting  the  load  impedance  of  the  receiving 
antenna. 

Figure  35  demonstrates  the  effect  of  the  relative  dielectric 
constant,  er,  on  the  transmitted  pulse.  An  increase  in  the  value  of 
the  relative  dielectric  constant  (er)  results  in  an  increase  in  magni- 
tude_of_the  received  waveform.  Since  the  characteristic  impedance 
*'uo/^oEr  is  decreased  by  an  increase  in  er,  an  increase  in  the  magnitude 
°f  the  received  voltage  is  to  be  expected.  Note  that  the  start  of 
the  waveform  for  tr  -  20  is  delayed  30  nsec  from  the  start  of  the  er  =  10 
waveform.  The  computer  program  is  also  accurately  predicting  the  delay 
time  of  the  transmitted  signal,  r  =  d/ygt  =  d/e r/c,  which  obviously 
increases  with  increasing  er.  Figure  35 b  displays  the  effect  of  er 
on  the  amplitude  spectra  of  the  transmitted  signals,  and  Fig.  35c  shows 
the  phase  spectra.  Note  that  the  amplitude  of  the  received  voltage 
for  er  =  20  is  lower  at  high  frequencies  (>50  MHz)  than  the  er  =  10 
case.  The  phase  is  also  delayed  by  increasing  er. 


a  = 
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Effect  of  dielectric  constant  on 
pulse  propagation. 

(a)  time  domain  waveforms 

(b)  amplitude  spectra 

(c)  phase  spectra 


The  previous  calculations  were  completed  assuming  a  bare  wire 
antenna,  which  is  really  not  the  case  when  the  antenna  is  resting  on 
the  absorber.  The  effect  of  an  insulating  dielectric  sleeve  on  the 
antenna  of  cr  =  1.0,  and  1/16"  thickness,  is  illustrated  in  Figs.  36 
and  37.  Figure  36  shows  the  transmitted  pulse  and  spectra  for  tr  =  10, 
a  =  .001  mho/m  and  d  =  20  feet  without  and  with  insulation,  and 
Fig.  37  shows  the  measurements  for  cr  =  10,  o  =  .01  mho/m  and  d  =  20  feet, 
without  and  with  ir.su.ation,  respectively.  Note  tne  differences 
between  the  insulated  and  uninsulated  cases  for  each  conductivity,  and 
compare  the  differences.  Obviously,  for  low  loss  media,  the  effect  of 
the  insulation  will  not  be  as  great  as  for  a  high  loss  medium.  Figures 
36  and  37  confirm  this  reasoning.  It  can  be  seen  from  the  o  =  .01  case 
(Fig.  37b)  that  the  coupling  between  the  bare  wire  antennas  is  stronger 
at  low  frequencies  while  the  coupling  for  the  insulated  antennas  is 
stronger  at  higher  frequencies. 

Figure  38  shows  that  for  constant  values  of  cr  and  a,  the  pulse 
shape  does  not  change  by  increasing  the  path  length.  Only  the  magnitude 
of  the  waveform  and  amplitude  spectrum  is  decreased.  There  is  also 
a  phase  delay  in  the  phase  spectrum,  which  corresponds  to  a  time  delay 
in  the  time  waveform. 

Another  important  factor  in  determining  the  received  waveform 
is  the  load  impedance  of  the  receiving  antenna.  To  illustrate  the 
dependence,  Fig.  39  compares  the  differences  between  pulses  transmitted 
on  the  same  antennas  in  the  same  medium  with  load  impedances  of  50q 
and  300q.  All  previous  waveforms  were  calculated  with  a  300U  load  to 
roughly  correspond  to  the  300q  twin  lead  of  the  measurement  system. 

The  impedance  of  the  twin  lead  is  really  lower  than  300^  when  it  is 
lying  on  the  interface  between  free  space  and  the  ground,  but  the  actual 
value  of  the  impedance  changes  with  both  frequency  and  the  parameters 
of  the  medium.  The  ideal  value  of  load  imoedance  would  naturally  be  a 
value  that  matches  the  impedance  of  the  antenna.  Note  from  Fig.  39 
that  a  change  in  a  purely  resistive  load  only  affects  the  magnitude 
of  the  received  waveform  and  the  amplitude  spectrum.  The  phase  spectrum 
remains  the  same. 

The  relative  dielectric  constant  and  conductivity  of  a  sub¬ 
surface  medium  are  not  constants  but  functions  of  frequency.  The 
relative  dielectric  constant  follows  a  general  decreasing  trend  with 
increasing  frequency,  while  the  conductivity  tends  to  increase  with 
increasing  frequency.  Just  what  changes  occur  in  the  values  is 
primarily  determined  by  the  moisture  content  of  the  medium.  Figure 
40  demonstrates  the  case  of  pulse  transmission  over  two  path  lengths 
in  a  dispersive  medium.  The  relative  dielectric  constant  and  con¬ 
ductivity  were  varied  linearly  over  the  2  to  40  MHz  ranqe  from  20  to 
10  and  .005  to  .015  mho/m,  respectively,  and  then  held  constant  at 
er  =  10  and  o  =  .015  from  40  MHz  to  100  MHz.  While  this  may  not  be  a 
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3&  r  dependence'  '*  is  helpful  simply  to  demonstrate  the 

evident  tL  'fl"6  “"Stitutive  parameters.  From  Fig.  40  it  is 
thv’d,nt  thab  the  base  width  of  the  time  domain  siqnals  is  areater  than 
the  base  width  of  the  constant  parameter  signals  (Fias  34^19 )  Thic 

rus  .rsss-.r*  -v r,!5iT5 

frequencies  accounts  for  the  spreading  effect  9 

nhtflinpH9^  compares  a  calculated  waveform  using  parameters  that  wer° 
o  tamed  from  the  propagation  measurements  throuqh  limestone  to  a 

SfaeUr^  10Uand  ar-nS0n5teh  1°  ^  °f  limestone-  Constant  values 

or  er  lu  and  o  -  .005  mho/m,  which  roughly  corresponded  to  the  HaSa  nf 

Fig.  14,  were  used  to  obtain  the  calculate  wavefbrE  When  caSoaSri  tJ 
re™rtedUu1?natr?hT!Ued  5U,S?  thrOU9t’  20  feet  of  "^tone,  wMch  was 

horizontal^'  'j"'6  IhoOT^orl^infllilesIml ' 

orizonta,  dipole  source  located  on  the  interface  that  the  horizontal 

cmly  iner  the  f^ds  P°  6’  h°We',er'  the  1nterface  does  5l9nif1- 

( Section  V)  nf^w*!  -=st  °f  accurac*  of  the  propagation  methods 

£  cEu&S  tS.lS2,(Siue(;)  f-s  pprd  to 

the^aluef  hJ0/m  resPfcctl've]y-  Fig.  42  shows  the  results.  rAboJe  iTmHz 
but  thc.UeS  2btrer  f0r  the  dielectric  constant  are  fairly  accurate 
,c  MJ  conductivity  is  slightly  inaccurate.  For  frequencies  below’ 

Ld  c  use  ethCpnpUCtin^  tdrm  °I  k2  =  begTns  to  domina?e 

inaccuracy  of  1211 1 3  T  ,As  the  ?onductivUy  is  increased,  the 
h„fCfh  aCyi°f  r  spreads  t0  frequencies  sliqntly  higher  than  15  MHz 
but  the  values  of  conductivity  becomes  more  accurate  ' 


VIII.  PROBE  CALIBRATION 

In  an  earlier  section  of  this  report  (Section  VIII  thp  flttpmutinn 
and  dispersion  effects  of  realistic  rock  media  on the  p/opagation  of 

sVtate  of  thp51^3  S  f re  demonstrated-  Calculations  were  made  using 

anafvs’il'^f  ^hi>ana  yS-S  3"d  Computer  Proqrems  for  the  complete 
ysis  of  arbitrary  wire  antennas  in  an  arbitrary  lossy  medium 

=■  n,rur  Mfes: 
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Fig.  40.  Effect  of  dispersion  on  pulse  propagation. 

(a)  time  domain  waveforms 

(b)  amplitude  spectra 

(c)  phase  spectra 
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be  ava’lable  from  the  National  Technical  Information  Service.  The 
program-,  have  been  made  available  for  our  use,  and  have  provided 
computational  advantages[l  ,2].  To  complete  a  theoretical  analvsis 
capability  for  the  pulse  sounding  probe,  a  correction  of  the  programs 
to  account  ^or  an  air-rock  interface  is  required.  This  correction, 
as  will  be  seen,  is  neither  obvious  nor  simple. 

Consider  first  what  is  meant  by  probe  calibration.  Basically  this 
means  that  the  transfer  function  of  the  dipole  antenna  on  the  surface  of 
the  medium  is  known.  Given  this  transfer  function  one  can  then  calcu¬ 
late  the  spectral  content  of  the  signal  actually  incident  on  the 
subsurface  target,  assuming  a  homogeneous  medium.  Note  that  this 
last  step  requires  a  knowledge  of  the  frequency-dependent  constitutive 
parameters  of  the  medium.  Knowing  the  two-way  effects  of  the  antennas 
and  the  medium,  the  response  waveform  spectrum  of  an  unknown  target 
can  then  be  normalized  to  obtain  the  spectrum  of  the  target's  impulse 
response  and,  via  fast  Fou-ier  transform  routines,  the  impulse 
response  of  the  target.  A  simplified  linear  system  analysis  of  the 
pulse  sounding  system  was  given  in  an  earlier  report[l].  The 
simplification  involved  was  that  the  transfer  function  of  the  trans¬ 
mitting  and  receiving  dipoles  were  assumed  to  be  identical  for  tne 
target  and  no  target  situations,  as  was  the  effects  of  the  medium. 

Of  course,  in  a  strict  sense  the  transfer  function  of  the  antennas 
cannot  depend  on  the  presence  or  absence  of  a  subsurface  target. 

However,  the  transfer  function  does  have  an  angular  dependence  and 
in  the  no  target  measurement  those  fields  which  penetrate  in  the 
direction  of  the  target  are  simply  not  obtained.  Similarly,  the  medium 
effects  will  be  somewhat  different  for  the  two  cases  because  different 
portions  of  the  medium  are  involved. 

Under  certain  conditions  which  will  not  often  bp  ^countered  in 
the  field,  it  is  possible  to  experimentally  estimate  the  transfer 
function  of  the  antenna  from  propagation  measurements .  The  tunnel 
geometry  in  limestone  described  earlier  provides  such  conditions.  Lump 
the  pulse  spectrum  and  that  of  the  cables,  baluns  and  twin  leads  as 
G(s),  then  if  two  identical  antennas  a^e  used  the  received  pulse 
spectrum  is  G(s)  Af(s)  M(s),  where  A-j  ( s )  is  the  transfer  function  of 
the  antenna  in  situ  on  the  rock  surface  and  M(s)  is  the  transfer  function 
of  the  medium.  From  earlier  propagation  measurements  over  two  different 
path  lengths  M(s)  is  known.  If  R (s )  is  the  received  pulse  spectrum 
then 

1/2 
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The  spectra  of  1  of  the  2  small  dipoles  shown  in  Fiq.  lc  are  shown  in 
Fig.  43.  These  data  were  obtained  using  the  Ikor  pulse  generator. 


A.  (s)  = 


R(s) 
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PT(s)  =  v G (s )  A,  (s ) 
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Note  that  the  currents  flowing  on  the  antenna  are  the  same  as  those 
when  the  antenna  is  immersed  in  an  infinite  lossy  medium  As  the 
antenna  is  brought  to  the  medium-air  interface,  the  image  approximation 
s  open  to  serious  questions.  However,  one  can  modify  the  currents 

“g  that  the  prop^gatUTctor 
ror  tnese  currents  is  that  for  an  infinite  wire  of  very  thin  radius 

ying  on  the  half-space.  Wait[12]  has  shown  that  the  first  term  of 

facto?sr°fo?aJi>n/HCJhr  !S  Simp1y  the  geometric  mean  of  the  propagation 
factorsfor  air  and  the  lossy  medium.  From  these  two  modifications, 

]:e-:  and  current,  estimates  of  the  antenna  input  impedance  and 

the  subsurface  fie  ds  can  be  obtained.  Calculations  of  this  type  are 
fields^the"?  W1  be.  1 compared  to  some  experimental  probing  of  the 
Ire  reaslnlble  96  P  determine  if  the  half-space  corrections 


IX.  CONCLUSIONS 

Significant  progress  ha:  been  made  in  the  experimental  develoDmpnt 
of  an  electromagnetic  pulse  sounding  probe  for  use  in  hard  rock  media 
A  full  scale  version  of  the  probe  has  been  built  and  its  effective  ’ 

Stailv  lnThproCt  r?dla,(1lmestOTe'  has  been  demonstrated  expert- 
mentally.  The  matching  does  not  require  disturbance  of  the  medium 

.e.,  the  ground  rods  used  for  earlier  overburden  tests[l]  have  been 
eliminated.  A  limiter  for  protection  of  the  sampling  head  used  as  a 
receiver  in  the  system  has  been  built  and  tested.  The  limiter  will 
permit  the  probe  to  oe  used  in  the  direct  mode  using  full  pulse  power. 

Using  elements  of  the  probe  and  2  different  pulse  generators 
propagation  measurements  through  up  to  28  feet  of  limestone  have  been 
.  .  *  Tom  data,  methods  for  obtaining  the  frequency-dependent 
attenuation,  relative  dielectric  constant  and  conductivity  of  the  rock 
medium  have  been  demonstrated.  From  these  same  data,  a  technique  for 

been '  "nus^ated  ,;empc,''a,  transfer  Actions  of  the  Srobe  has 

n  .h  Ful]  sc5le  fc'Perimental  measurements,  using  the  probe  in  an 
orthogonai  mode,  have  been  rede  of  fault,  joint  and  lithologic  con- 

45  flet  ?oTlSl?thCnim-9  ln  J  d°l0niite  qt'arry'  The  raximum  range 
45  feet  lor  a  lithologic  contrast  target.  This  is  not  necessarily  the 

maximum  range  of  the  electromagnetic  sounder,  however,  -  present 

In^fvtl0nS  alVhat  this  range  is  somewhat  ^ss  than  80  feet.  As  with 
any  experimental  measurement  of  an  uncontrolled  and  inaccessible  tarqet 
the  interpretation  given  the  assured  data  is  always  subject  to 
qualification.  It  is  maintained,  however,  that  internal  reflection 
Sni5HS’  mUltip1e  reflectio"s  on  the  cables  etc.,  have  been 

?  e  IeasLand  sUe'I  F  n?!SUred  ^  f1t  the  ^logical  descriptions  of 
slt®s\  Furthermore,  none  of  the  measured  data  are  incon¬ 
sistent  with  geological  interpretations. 
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.  Theoretical  advances  have  been  modest  because  of  the  emphasis 
during  toe  present  contract  period  on  experimental  measurements 
Realistic,  pulse  propagation  calculations  employing  finite  antenna 
elements  iri  a  lossy  medium  with  frequency  dependent  constitutive 
parameters  have  been  made.  These  are  not  quasi -static  or  idealized 
dielectric  approximations,  but  nurren'cal  calculations  where  the  sub¬ 
stitution  of  assumed  or  measured  frequency-dependent  constitutive 
parameter  for  the  medium  are  the  only  real  approximation.  Some 
progress  has  been  made  to  extend  these  calculations  to  include  an  air- 
medm  interface  but  the  analyses  are  incomplete  as  yet.  The  problem 
is  compounded  by  the  fact  that  for  realistic  rock  media,  the  media  is 
both  a  conductor  and  dielectric  for  the  video-type  pulses  of  interest. 
Pulse  propagation  calculations  which  include  an  air-medium  correction 
will  be  made  during  the  present  contract  period.  The  validity  of  the 
air-medium  correction,  however,  may  not  be  fully  explored 


.  is1felt  that  at  this  tlme  the  electromagnetic  pulse  sounding 
probe  is  fully  ready  for  testing  in  an  actual  hard  rock  tunnel  site. 

The  system  at  present  is  a  laboratory  instrument,  i.e.,  components  are 
not  necessarily  field-hardened,  but  permit  measurements  in  the  field 
Successful  measurements  reported  here  fully  justify  testing  of  the  pulse 
sounding  probe  in  an  actual  tunnel  geometry.  A  fully  calibrated  probe 
and  system  would  not  be  necessary  at  present.  Formal  testing  of  the 
electromagnetic  pulse  sounding  probe  is  recommended  to  provide  a 
realistic  comparison  with  current  acoustic  and  seismic  techniques. 


X.  FUTURE  PLANS 

During  the  second  interim  of  the  contract  it  is  intended  to 
concentrate  on  2  primary  objectives:  Refinement  of  the  antenna  analysis 
to  permit  half-space  corrections  to  be  made  for  an  arbitrarily  shaped 
wire  antenna,  and  a  second  measurement  program  at  the  Plum  Run  Quarry 
site.  The  new  measurements  will  include  direct  as  well  as  orthogonal 
mode  responses  using  peak  pulse  outputs.  In  addition,  the  original 
intention  of  interrogating  joint  and  fault  targets  with  the  probe  on 
vertical  rock  faces  will  be  completed.  This  will  permit  the  joint 
and  fault  targets  to  be  viewed  with  orientations  with  respect  to  the 
probe  which  are  more  likely  to  be  encountered  in  an  actual  tunr°linq 
operation. 

The  refinement  of  the  antenna  analysis  will  permit  calibration 
ot  the  probe  and  via  this  calibration  an  improved  processing  of  the 
target  responses.  In  addition,  a  design  capability  for  theoretically 
matching  the  antenna  geometry  to  an  arbitrary  medium  will  be  obtained. 
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